o 

O 

Q 



o 

s 

-I— > 

a 
a 

o 
o 



> 

O 
O 

<N 



1 



mito calculation of the local magnetic moment in titanium doped zinc 
oxide with a corrected-band-gap scheme 
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The local magnetic moment of Ti:ZnO is calculated from first principles by using the corrected-band-gap 
scheme (CBGS). The result shows that the system is magnetic with the magnetization of 0.699 /ie per 
dopant. The origin of the local magnetic moment is considered to be the impurity band partially occupied 
by the donor electrons in the conduction band. Further, the impacts of applying Hubbard U to Ti-d orbital 
on the magnetic moment have been investigated. 
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Zinc oxide, a diamagnetic semiconductor, doped by 
certain transition metals (TM) was observed to be a 
room- temperature ferromagnet. Moreover, according to 
the experiment of Venkatesan et a/., 1 "there is a strik- 
ing systematic variation (of magnetization), with max- 
ima near the beginning of the 3d series at Ti and V (0.5 
Mb)-". 

To understand the origin of the magnetic moment, 
first-principle calculations for TM doped ZnO have been 
carried out, especially for cobalt. However, there are 
few reports on Ti doped ZnO. Early calculation 2 us- 
ing density functional theory (DFT) approach based 
on GGA (Generalized Gradient Approximation) and 
GGA+U (GGA plus on-site Coulomb repulsion) pre- 
dicted a magnetic moment (0.9 /ie and 1.6 /ie, respec- 
tively) larger than the experimental result 1 . In addition, 
structure optimization was ignored in Ref. 2. Further 
work 3 with atomic-position relaxation using full poten- 
tial linear augmented plane wave method obtains a result 
of 0.63 /iB per dopant, which is in good accordance with 
the experiment. 1 However, both of the previous works 
may be affected by the well-known problem of local-spin- 
density approximation (LSDA)/GGA, underestimation 
of the band gap. As a result, they are suffered from the 
incorrect impurity levels relative to the host band edges, 
which may significantly impact the calcualted local mag- 
netic moment originating from the impurity levels. In 
this work, we investigate Ti:ZnO using a corrected-band- 
gap scheme (CBGS), where both band gap and impurity 
levels are corrected by the DFT+ U approach 4 . 

A 3 x 3 x 2 supercell with 72 atoms is constructed 
from the wurtzite primitive cell of ZnO with the exper- 
imental lattice parameters, and one Zn atom is substi- 
tuted by Ti. The atomic positions are first optimized by 
using Perdew-Burke-Ernzerhof (PBE) 5 parameterization 
of GGA and projected augmented wave (PAW) method. 
The optimization is accomplished by using VASP code 6 



TABLE I. The lengths (A) of Ti-0 bonds. The unique refers 
to the bond along the c axis, and the other refers to those 
bonds in the basal plane. 





unique 


other 


initial 


1.99 


1.97 


optimization 


1.88 


1.88 



on a 3 x 3 x 3 A:-mesh. It shows that the four O atoms 
nearest neighbor to Ti move inward. As shown in Table HI 
the lengths of the four Ti-0 bonds, one along the crys- 
tallographic c axis (the unique) and others in the basal 
plane, are all shortened to 1.88 A, which are originally 
1.99 A and 1.97 A in the primitive cell of ZnO, respec- 
tively. As a result, the local symmetry changes from Cs v 
toT d . 

After structural optimization, the electronic structure 
and magnetic properties of Ti:ZnO are calculated with 
the CBGS, by which the Hubbard Us applied to Zn-d 
orbital (5 eV) and O-s orbital (-50 eV) enhance the band 
gap of pristine ZnO to 2.8 eV 7 . The A:-mesh increases to 
4x4x4 and the plane- wave cutoff is 500 eV. Furthermore, 
with CBGS, calculation with Hubbard U applied on Ti-d 
orbitals (UTi-d) ar e performed to investigate the impacts 
of ^-electron correlation. 

In Table [III we summarize the results of magnetic mo- 
ments by different methods. Obviously, the magnetic 
moments are increased by CBGS and further enhanced 
by CBGS+ U, and reach the largest value when U =3 eV. 
The total density of states (DOS) [Fig. Ufa)] shows that 
the gap between the O-p band and Ti-d band is pre- 
served in Ti:ZnO with CBGS. The valence band (VB) of 
the host and the Zn-d band are pushed down to lower en- 
ergy by CBGS. Moreover, the hybrid between Zn-d band 
and O-p band is reduced. These changes due to CBGS 
are observable in Co: ZnO as well 7 . As shown in Figs. [1] 
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FIG. 1. (a) The total DOS by GGA and CBGS, and the 
projected DOS of Ti-3d (b) by GGA and (c) by CBGS. 
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FIG. 2. Formation of defect levels in TiiZnO with (, 
and (b) CBGS+£/ T i-d. The Ti|+ orbitals in tetrahedral site 
splits into two- fold e and three- fold £ levels in crystal field. 
Then the t levels interact with the dangling-bond states of 
the Zn vacancy, and form the defect levels. 



TABLE II. The magnetic moment given by different methods. 



Magnetic moment (/xb) 



Method 


Total 


Ti 


U T i-d (eV) 


GGA 


0.227 


0.199 




CBGS 


0.699 


0.703 




CBGS+U T i-d 


0.777 


0.809 


2.0 




0.973 


0.957 


3.0 



(b) and (c), CBGS increases the exchange split of Ti-d 
bands from 0.25 eV to 0.75 eV. The majority Ti-d band 
lies at the conduction-band minimum (CBM) and strad- 
dles the Fermi level. Thus, the system is metallic. 

To understand the origin of the local magnetic mo- 
ment, we introduce a model about the formation of the 
defect levels (Fig. [2]). According to Zunger et al^ the 
defect levels are a result of the hybridization between the 
levels of the host cation vacancy (anion dangling bonds) 
and the 3d orbitals of the impurity under the ligand field. 
As shown in Fig. [2j we provide the diagram for Ti sub- 
stitution in ZnO. Because of the conductive character of 
Ti-d electrons, the defect levels, two hybrid t levels and 
two nonbonding e levels, lie above CBM. Since the e+ FR 
is occupied by 2e, the magnetic moment should be 2 /iB 
by the first Hund's rule. Fig. [2] also illustrates the im- 
pacts of adding U on Ti-d on the defect levels, that the 
exchange split energy is enhanced by applying U, and 
that the minority defect level e^ FR moves to the energy 
higher than the majority defect level £^ FR . 

According to the above analysis (Fig. [2]), the total mag- 
netic moment should be 2 /iB per dopant. However, these 
values in Table HJ are smaller than 1 /iB- The deviation is 
caused by the band filling effect. As shown in Fig. [3ja), 
the defect levels form the impurity bands. As the low- 
est impurity band occupied by two Ti-d electrons (donor 




(a) (b) 

FIG. 3. The impacts of band filling effect on magnetic mo- 
ment. The occupied impurity band is higher than the host 
CBM. Consequently, part of the donor electrons transfers 
from the impurity band to the host CBs, and the Fermi energy 
lowers down. 



electrons) are above the host CBM, part of the donor 
electrons transfers from the impurity band to the host 
CBs, and the Fermi energy lowers down [Fig. [3fb)] . In 
addition, we integrate the DOS of CBGS (Fig. [I]) from 
CBM to Fermi level and obtain the electronic popula- 
tions of 1.4 and 0.6 for the majority and minority spins, 
respectively. Hence, the total number of donor electron 
is 2.0, supporting the defect level model (Fig. [2} and the 
band filling model (Fig. [3]). 

By applying Hubbard U on Ti-3d orbital, the exchange 
split of Ti-d band continues to increase and the positions 
of the impurity bands in different spins support the model 
in Fig. E(b), as shown in the DOS with CBGS+Z7 Ti -d 
(Figs. 3] and [5]). When U =3.0 eV, the double degenerated 
and the partially occupied impurity band corresponding 
to the e+ FR level splits into two bands, one of which is 
occupied by le and results in the magnetic moment of 
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FIG. 4. The projected DOS with CBGS+£/ T i-d (2.0 eV). 
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FIG. 5. The projected DOS with CBGS+*7 T i-d (3.0 eV). 



about 1 fiB (Table HI)). In general, a well-accepted value 
of U for Ti-3d orbitals is 2.0 eV— . At the same time, we 
also observe the changing of the spin density associated 
with the split of degenerated e orbitals as U increases 
from 2 to 3 eV (Fig.©. When {7=2.0 eV [FigEJa)], the 
shape of spin density is obviously a combination of two 
d orbitals, d 3z 2_ r 2 and d x i_ y i. However, that of {7=3.0 
eV [FigEJb)] is just one of the two d orbitals, d x 2_ y 2. 

In conclusion, we calculate the local magnetic moment 
of Ti:ZnO from first principles with CBGS. The result 
shows that the system is magnetic with 0.699 {ib per 
dopant. The origin of the local magnetic moment is also 
discussed. We conclude that, due to the position of the 
defect levels and band filling effect, the Ti-d band is par- 
tially occupied by donor electrons. Further, we inves- 



tigate the impacts of adding U to Ti-d orbitals on the 
magnetic moment of the system. Both DOS and spin 
density map indicate that increasing {7xi-d leads to a 
further split of Ti-d orbitals and simultaneously causes a 




FIG. 6. The spin density maps from two respects with 
CBGS+{7 T i-d, (a) {7=2.0 eV and (b) {7=3.0 eV. When {7=2.0 
eV, the shape of spin density is obversely a combination of two 
d orbitals, d 3z 2_ r 2 and d x 2_ y 2. However, that of {7=3.0 eV 
is just one of the two d orbitals, d x 2_ y 2. 



larger magnetic moment. 
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